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ABSTRACT 

We report on an analysis of XMM-Newton data from the neutron star low 
mass X-ray binary (LMXB) Serpens X-l (Ser X-l). Spectral analysis of EPIC 
PN data indicates that the previously known broad iron Ka emission line from 
this source has a significantly skewed structure with a moderately extended red 
wing. The asymmetric shape of the line is well described with the laor and 
diskline models in XSPEC and strongly supports an inner accretion disk origin 
of the line. To our knowledge this is the first strong evidence of a relativistic line 
in a neutron star LMXB. This finding suggests that the broad lines seen in other 
neutron star LMXBs likely originate from the inner disk as well. Detailed study of 
such lines opens up a new way to probe neutron star parameters and their strong 
gravitational fields. The red wing of the iron line from Ser X-l is not as broad 
as that observed from some black hole systems. This is not unreasonable for a 
neutron star system, as the accretion disk has to terminate at or before the hard 
stellar surface. Finally, the inferred source inclination angle in the approximate 
range 40°-60° is consistent with the lack of dips and eclipses from Ser X-l. 

Subject headings: accretion, accretion disks — line: profiles — relativity — stars: 
neutron — X-rays: binaries — X-rays: individual (Serpens X-l) 



1. Introduction 

The bright persistent LMXB Ser X-l was discovered in 1965 (Bowyer et al. 1965). The 
detection of type I X-ray bursts from this source established that it harbors a neutron star 
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(Swank et al. 1976; Li et al. 1977), as such bursts originate from the thermonuclear burning 
of matter on the stellar surfaces (Woosley & Taam 1976; Lamb & Lamb 1978; Strohmayer & 
Bildsten 2006). The lack of energy dependent dips and eclipses in the X-ray light curve of Ser 
X-l suggests that its inclination angle may be less than 60° (Frank et al 1987; White & Swank 
1982). The analysis of RXTE and BeppoSAX data, revealed that the continuum spectra could 
be well fitted with a combination of absorbed Comptonization and disk blackbody models 
(Oosterbroek et al. 2001). Moreover, these authors reported the existence of a broad iron 
emission line near 6 keV, which they could adequately fit with a Gaussian model. Broad 
iron emission lines (near 6 keV) have been observed from many LMXBs (Asai et al. 2000; 
Bhattacharyya et al. 2006; Miller et al. 2002a), and active galactic nuclei (AGN; Reynolds & 
Nowak 2003). These lines could be broadened by either Doppler and relativistic effects (due 
to Keplerian motion in the inner accretion disk; Fabian et al. 1989), or Compton scattering 
(in a disk corona; Misra & Kembhavi 1998). It has also been suggested recently that complex 
absorption from perhaps several ionized components may partially account for some of the 
broad line component (e.g., Turner et al. 2005). The ASCA data from the Seyfert-1 galaxy 
MCG-6-30-15 gave the first strong clues that such lines originate in the inner accretion disk, 
as the high signal-to-noise ratio data revealed a skewed and double-horned profile, consistent 
with Doppler and relativistic effects in the inner accretion disk, where the speed of matter 
is a substantial fraction of the speed of light (Tanaka et al. 1995). Analysis of more recent 
data has strongly supported this inner disk origin of the line for MCG-6-30-15 and other 
AGN (Reynolds & Nowak 2003; Wilms et al. 2001; Reynolds et al. 2004; Reeves et al. 2007; 
Miniutti et al. 2007). Such lines, therefore, are extremely useful probes of strong gravity, 
and black hole properties, including spin (Brenneman & Reynolds 2006). 

An inner disk origin has also been suggested for broad lines in some stellar mass black 
hole binaries, such as Cyg X-l and GRS 1650-500 (Miller et al. 2002a; 2002b; Fabian et 
al. 1989). But, the origin of the broad iron lines from neutron star LMXBs has not been 
as well understood. This is largely because the available data have a modest signal-to- 
noise ratio, which was not precise enough to rule out simple symmetric profiles (such as a 
Gaussian profile). If the inner accretion disk origin can be established for broad iron lines 
from neutron star LMXBs, then the shape and width of the line may be used to constrain 
the radius of the inner edge of the disk, as well as the Keplerian speed at that radius. The 
former can give an upper limit to the neutron star radius (as the disk inner edge radius 
must be greater than or equal to the stellar radius), while both quantities may be utilized to 
constrain the stellar mass. Moreover, high frequency quasi-periodic oscillations (kHz QPOs) 
are observed from many neutron star LMXBs (van der Klis 2006). Although the actual 
physical mechanism responsible for these timing features is still debated, the leading models 
involve strong gravity, with the accretion disk close to the neutron star surface. Therefore, 
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if the inner disk origin can be established, the broad iron lines can be useful to constrain 
kHz QPO models, and together these two spectral and timing features will be extremely 
important as probes of strong gravity, and to constrain neutron star parameters. We note 
that the constraints on neutron star mass and radius can be useful to determine the equation 
of state (EoS) of high density cold matter at the stellar core, which is a fundamental problem 
of physics (see Bhattacharyya et al. 2005). In this Letter, we present evidence for an inner 
accretion disk origin of the broad iron line from the neutron star LMXB Ser X-l. In § 2, we 
describe our analysis of XMM-Newton data, and in § 3 we discuss the implications of our 
results. 



2. Spectral Analysis 

The neutron star LMXB Ser X-l was observed three times with XMM-Newton in March, 
2004 (obslds: 0084020401, 0084020501, 0084020601). These observations were each sepa- 
rated in time by two days, and each consists of about 22 ks of data. Due to the brightness of 
the source, there are no European Photon Imaging Camera (EPIC) MOS detector data avail- 
able for these observations. EPIC PN data were obtained in timing mode, and the source 
photons were also registered by the Reflection Grating Spectrometer (RGS) instruments. 
However, as no significant source spectral lines were detected in the RGS energy range, here 
we report only the analysis of the EPIC PN spectra. We have excluded the portions of the 
EPIC PN data with high and variable background (due to soft proton flares), and extracted 
source spectra, background spectra, and response matrices using the XMM-Newton Science 
Analysis Software (SAS; version 3.0). The SAS task 'epatplot' has indicated that this set of 
spectra (set 1) is modestly piled up. We have, therefore, extracted another set of spectra (set 
2) after excluding the contributions from the central (and hence the brightest) pixels in order 
to minimize the pile-up effect. Results from the two spectral sets are consistent with each 
other, and the broad double-peaked line, which is the focus of this Letter, is present in both 
sets in the same energy range. This is expected, as pile-up cannot generate a double-peaked 
line. Based on these findings we do not believe pile-up is a significant concern for our study. 
Therefore, we have used spectral set 1 because of its higher signal-to-noise ratio. 

We have rebinned the spectra (see, for example, Ibarra et al. 2007), fitted them with 
various models (within XSPEC), and found that the best model to describe the observed 
continuum spectra consists of an absorbed Comptonization (compTT) component plus a disk 
blackbody (diskbb). The addition of diskbb to the compTT component is essential, as can 
be seen from Table 1. This is in accordance with observations of Ser X-l obtained with 
other X-ray missions (see Oosterbroek et al. 2001). Oosterbroek et al. noted that the choice 
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between a disk blackbody and a single temperature blackbody was arbitrary, however, our 
spectral fitting shows that a disk blackbody is preferred to a simple blackbody. For example, 
the choice of the former reduces x 2 by about 47 compared to that for the latter (obs. 1, i.e., 
obsld 0084020401). Such a decrease in \ 2 is significant, and a similar decrease is found for 
obs. 2 (obsld 0084020501) and obs. 3 (obsld 0084020601). Table 1 shows that the addition 
of a Gaussian emission line near 0.54 keV improves the fit very significantly. As we have 
not found any source spectral lines at this energy in the RGS data, we suspect that this 
feature may be of instrumental origin, so we have fixed the line centroid of this feature at 
~ 0.54 keV in our subsequent model fitting. Finally, we have found highly significant excess 
emission near 6 keV. As a broad iron emission line is detected near this energy from many 
LMXBs (see § 1), we have fitted this excess emission with a Gaussian. This extra component 
is very significant (see Table 1), consistent with the results reported in Oosterbroek et al. 
(2001). However, as we have mentioned in § 1, it is not yet known with any certainty 
what broadens this line for neutron star LMXBs. To explore this question in more detail 
we have replaced the 6 keV Gaussian component with the diskline model component of 
XSPEC. This component represents the spectral line emission from the inner portion of an 
accretion disk in the Schwarzschild spacetime (Fabian et al. 1989). The corresponding fits 
show that the diskline model describes the line profile significantly better than the simple 
Gaussian profile (see Table 1). This result supports the notion that the line is produced in 
the inner accretion disk. In order to further substantiate this result, we have also modeled 
the spectrum using the laor model. This model is similar to the diskline model, but it 
includes the effects of the spin of the central star (Laor 1991). The fact that the laor model 
also fits the line profile significantly better than the simple Gaussian (see Table 1) provides 
strong evidence of the inner disk origin of this line. In fact, laor describes the line profile 
slightly better than diskline (primarily for obs. 1; see table 1). However, while there is 
some indication of a smaller disk inner radius in the laor fit using obs. 1, a value of 6r g 
is not yet excluded with sufficient significance to argue that the spin of the neutron star is 
affecting the line profile. In principle, deeper observations could test for such an effect. 

In Table 2, we give the best-fit parameter values for the XSPEC model wabs* (compTT+ 
diskbb+gauss+laor) for each observation. These values are consistent across the three ob- 
servations. Here we note that we have used the disk geometry for the Comptonizing plasma 
(compTT model). However, even the spherical geometry gives very similar parameter values, 
except for the compTT optical depth tq (the latter geometry gives tq ~ 15). We have con- 
strained the rest frame line energy (E^, laor component) in the range 6.4 — 6.97 keV while 
fitting. This is because, Ka spectral lines from neutral or ionized iron are expected in this 
range. However, although we have always found 6.4 keV as the best-fit value of E L , relaxing 
this constraint does not lower the best-fit E L much. We have found the best-fit source inch- 
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nation angle, zl, in the range as 40° — 60°, which is consistent with the expected value (< 60°) 
for Ser X-l (due to the lack of eclipses and dips; see § 1). Here we note that, since there 
are many model parameters, we have fixed the line parameters (except the normalizations) 
to their best-fit values to estimate the errors on the parameters of the wabs, compTT, and 
diskbb components. Similarly, we have frozen the parameters (except the normalizations) 
of the compTT and diskbb models in order to estimate the errors on the spectral line param- 
eters. Therefore, the quoted error values are likely somewhat underestimated. Fig. 1 shows 
the data, model components, and data-to-model ratio for observation 2. The laor model 
component for the broad iron line is clearly shown by the double-peaked dotted line. The 
moderately structured shape of the data-to-model ratio results in a relatively high overall 
reduced x 2 f° r Ser X-l (Table 2). Note that similar residuals, and a high reduced x 2 were 
also reported by Oosterbroek et al. (2001; see their Figs. 5 & 6). These authors argued 
that these narrow structures could not be caused by an incorrect continuum modeling, since 
such incorrectness would probably give more smoothly varying data-to-model ratios. Fig. 2 
exhibits the structure of the broad relativistic iron line. The two panels show this line from 
two observations. The data points of the figure clearly show the line with an extended red 
wing. This figure also explicitly shows that both the laor model and the diskline model 
(dotted line in each panel) fit the spectral line profile well. 



3. Discussion and Conclusions 

In this Letter, we report the results of the spectral fitting of the XMM-Newton EPIC 
PN data from the neutron star LMXB Ser X-l. The best-fit continuum parameter values 
(see Table 2) are generally consistent with those of Oosterbroek et al. (2001). Our slightly 
higher Comptonization plasma temperature (T c ) compared to that of Oosterbroek et al. can 
be explained in terms of the lower 2 — 10 keV flux that we have found. This is because, as 
the intensity of an LMXB decreases, its energy spectrum generally becomes harder. 

It was previously known that Ser X-l exhibits a broad iron line (see § 1), but earlier 
data did not have the statistical quality to rule out simple symmetric profiles (such as a 
Gaussian). As a result, its origin was not strongly constrained. Knowledge of the line's 
origin is very important, because if it is produced in the inner accretion disk, it can be used 
(1) to probe the strong gravitational field of the neutron star, (2) to constrain the stellar 
parameters, and (3) to constrain models of the kHz QPOs, which may, in turn, be useful to 
achieve the first two goals (see § 1). By analyzing XMM-Newton EPIC PN data we have 
demonstrated that the line is significantly asymmetric, and can be accurately modeled with 
physical models of line formation in the inner disk (using the diskline and laor models 
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in XSPEC; see Fig. 2). This strongly supports the idea that this spectral line originates 
from the inner accretion disk of Ser X-l (see § 1). The high signal-to-noise ratio data also 
revealed the detailed shape of the broad and relativistically skewed line with an extended 
red wing (Fig. 2). Such lines have so far been observed from AGNs and a few Galactic black 
hole X-ray binaries (Tanaka et al. 1995; Miniutti et al. 2007; Miller et al. 2002a; 2002b), 
but this is the first strong evidence of a relativistic line in a neutron star system. As we 
have mentioned in § 1, other neutron star LMXBs exhibit broad iron emission lines, which 
so far have been adequately modeled with simple Gaussian profiles. Our finding for Ser X-l 
suggests that the lines seen in other LMXBs may also originate from the inner accretion disk, 
and its relativistically skewed nature could be confirmed with sufficiently long observations. 
If future observations bear this out it will open up an exciting new opportunity to probe 
neutron stars and strong gravity with deep spectroscopic observations. 

We have found that the equivalent width (EWl) of the iron line from Ser X-l is lower 
than that reported in Oosterbroek et al. (2001), which may suggest that the strength of 
this line decreases with the source intensity. The source inclination angle inferred from both 
the diskline and laor fits is consistent with < 60°, which is expected for the non-dipping 
nature of Ser X-l (see § 1). The best-fit values of R m for obs. 1 & 3 (see Table 2) suggest 
that the accretion disk extends almost to the neutron star surface. Although this is not the 
case for obs. 2, it is to be noted that both the source count rate and the iron line strength are 
relatively low for this observation (Table 2), and hence the spectral line statistics are not as 
good as for the other two obsIDs. However, we cannot make a definite conclusion regarding 
the disk extension with the current data. But we note that the accretion disk can plausibly 
approach the neutron star surface because the stellar magnetic field in LMXBs is relatively 
low (^ 10 8 ~ 9 G), and the implied accretion rate for Ser X-l is high (^ 0.28 — 0.32MEdd; 
based on the observed flux, and the Eddington luminosity ~ 2.0 — 3.8 x 10 38 ergs s -1 , for a 
distance of ~ 9.5 — 12.7 kpc; Jonker & Nelemans 2004). An important difference between 
neutron star and black hole systems in the context of disk line formation is the presence of 
the neutron star surface. This sets a firm limit on the inner edge of the disk in a neutron 
star system, even a rapidly spinning one, although the exact limit will depend on the EoS 
of neutron star matter. For black holes, significant spin always implies a smaller disk inner 
edge radius (for a corotating disk). Therefore, for very fast spinning black holes, the red 
wing of the relativistic iron line is expected to be very broad. So the relatively modest 
extended red wing (compared to some black hole systems; Miniutti et al. 2007) of the iron 
line from Ser X-l appears consistent with the fact that this source harbors a neutron star. 
This relativistic iron line is expected to be accompanied by a continuum reflection spectral 
component (Reynolds & Nowak 2003). A distinct feature of this component is a broad hump 
near 30 keV. As the signal-to-noise ratio of the EPIC PN data is small at higher energies 
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(> 10 keV), we could not determine if this component is present in the Ser X-l spectra. 

The authors thank Tim Kallman and Jean Cottam for useful discussions, and John 
Miller and an anonymous referee for helpful comments. 
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Table 1. Fitting of the XMM-Newton EPIC PN energy spectra from Ser X-l with various 



XSPEC models. 


No. 


XSPEC Model 


X 2 /dof a 
(Obs. l b ) 


X 2 /dof a 
(Obs. 2 b ) 


X 2 /dof a 
(Obs. 3 b ) 


1 


wabs*compTT 


2689.1 
706 


2517.2 
706 


2464.3 
706 


2 


wabs* (compTT+diskbb) 


1708.9 (Q) 
704 v u / 


^(0) 




3 


wabs* (compTT+diskbb+gauss) 


(8.0E-19) 


n 7 ° 3 2 8 (3.3E-13) 


13 7 3 8 2 2 (4.2E-11) 


4 


wabs* (compTT+diskbb+gauss+gauss) 


1J §T ( 2 - 7E - 36 ) 


^ (9.7E-23) 


10 6 9 9 8 9 ° (8.4E-30) 


5 


wabs* (compTT+diskbb+gauss+diskline) 




W (3-0E-5) 


W (3-1E-4) 


6 


wabs* (compTT+diskbb+gauss+laor) 


^ (3.6E-5) 


^ (3-4E-5) 


W (2-3E-5) 



a The number in the parentheses is the probability (calculated from F-test using XSPEC) of the 
decrease of x 2 /dof by chance from the value of the previous row to that of the current row (for models 
2 — 5). For model 6, the comparison is with model 4. A very small value is replaced with a zero. 

b Obs. 1: Obsld 0084020401; Obs. 2: Obsld 0084020501; Obs. 3: Obsld 0084020601. 
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Table 2. Best fit parameters (with 90% confidence) for the XMM-Newton EPfC PN 

energy spectra from Ser X-l. 



Parameter* 


Obs. l b 


Obs. 2 b 


Obs. 3 b 


iV H c (10 22 atoms cm" 2 ) 
T d (keV) 
T c c (keV) 

T f 

T in s (keV) 
E L h (keV) 

RiJ 
-Rout k 

(degree) 


_ . _ i n no 

0.48 + °;°2 

„ —_i_n 01 

0-1718:81 

q i o+u.04 
°- lo -0.36 
a I q+0.68 

D.iy_ 074 

i 71+O.IO 
' ± -0.12 

6.4oi8:8 s 

9 I 7+0.30 
-0.13 

4 04+ 2 - 14 

^• u ^-0.68 
971 c;+il5.2 
z ' i - d -93.7 

44.71^ 


_ . . 1 n ni 

u -^-o.oi 

o.i9l8:8i 

_ _ _ in 1 a 

z - uo -o.n 

7 9C+0.30 
' • z " J -0.32 
1 91+O.O4 
± - zi -0.04 

« i n+0 03 

6.4018:8' 

9 qO+0.31 

z - oo -0.56 

le.iol 8 ^ 

400.0l° 70 ° 9 
50.2±|;t 


_ . _ in no 

u -^ y -0.02 

0.1718:8? 

^• 8i -0.18 
a 71+O.48 
- 1 - 0.52 

1 4n+o-08 

i -^ u -0.08 
« 1 n+n 04 

6.4oi8:8 4 

9 rc+0.39 
z - Jd -0.30 

J - jy -1.19 

qi q r+86.5 
0±0.d_ 166 g 

39.71^ 


EWi m (eV) 


105.1118 


85.9l? 9 8 


95.2™ 


xVdof 


1153.6/696 


917.8/696 


1060.5/696 


Flux 11 (0.5 - 2 keV) 


0.93 


0.87 


0.91 


Flux 11 (2 - 10 keV) 


4.22 


3.34 


3.79 


Parameters of 


the 


XSPEC 


model 



wabs* (compTT+diskbb+gauss+laor) . 

b Obs. 1: Obsld 0084020401; Obs. 2: Obsld 0084020501; Obs. 
3: Obsld 0084020601. 

c Hydrogen column density from the wabs model component. 

d Input soft photon (Wien) temperature of the compTT model 
component. 

c Temperature of the Comptonizing plasma (compTT model 
component; disk geometry). 

f Optical depth of the Comptonizing plasma (compTT model 
component). 
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g Temperature at inner disk radius from diskbb model compo- 
nent. 

h Rest frame energy of the broad relativistic iron emission line 
(laor model component). 

'Power law index of emissivity (laor model component). 

j Inner radius (in the unit of GM/c 2 ; M is the mass of the 
neutron star) from the laor model component. 

k Outer radius (in the unit of GM/c 2 ) from the laor model 
component. 

'Source inclination angle from the laor model component. 

m Equivalent width of the broad relativistic iron emission line 
(laor model component). 

"Observed flux in 10 -9 ergs cm" 2 s -1 . 
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Energy (keV) 



Fig. 1. — XMM-Newton EPIC PN energy spectrum from Ser X-l for obs. 2 (obsld 
0084020501). The upper panel shows the data (error bars), model (solid line), and indi- 
vidual additive model components (dotted lines). Here we have used the best-fit XSPEC 
model wabs*(compTT+diskbb+gauss+laor) (see Table 2). Among the two continuum ad- 
ditive model components, the upper dotted line shows compTT, and the lower dotted line 
shows diskbb. The low energy emission line is the gauss component, while the high energy 
broad iron emission line is the laor component. The lower panel shows the data to model 
ratio. 
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Fig. 2. — Upper panel: XMM-Newton EPIC PN energy spectrum from Ser X-l for obs. 1 
(obsld 0084020401). With the error bars we show the data intensity in excess of the model 
intensity, and normalized by the model intensity (i.e., (Data— Model) /Model). Here the 
model is the best-fit model (of Table 2) minus the laor component. The laor component 
is separately shown with the dotted line. Lower panel: Similar to the upper panel, but for 
obs. 3 (obsld 0084020601), and for the diskline component (model 5 of Table 1) instead of 
the laor component. For each panel, the data points clearly show a broad relativistic iron 
emission line with an extended red wing. 



